Dynamics of a quantum quench in an ultra-cold atomic BCS superfluid 
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We study dynamics of an ultra-cold atomic BCS superfluid driven towards the BCS superfluid-Fermi liquid 
quantum critical point by a gradual decrease of the pairing interaction. We analyze how the BCS superfluid falls 
out of equilibrium and show that the non-equilibrium gap and Cooper pair size reflect critical properties of the 
transition. We observe three stages of evolution: adiabatic where the Cooper pair size is inversely proportional to 
the equilibrium gap, weakly non-equilibrium where it is inversely proportional to the non-equilibrium gap, and 
strongly non-equilibrium where it decouples from both equilibrium and non-equilibrium gap. These phenomena 
should stimulate future experimental characterization of non-equilibrium ultra-cold atomic BCS superfluids. 
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Ultra-cold atoms provide quantum simulators for emulat- 
ing challenging theoretical models. This offers invaluable op- 
portunities for the study of many-body systems which used 
to be experimentally inaccessible by conventional condensed 
matter setups 1 1 ] . The ability to make the parameters of the 
emulated models time-dependent encourages exciting stud- 
ies of non-equilibrium dynamics of many-body quantum sys- 
tems This opens up fresh prospects for interdisciplinary 
research of dynamics of phase transitions, where universal- 
ity of the critical behavior links dynamics of distinct systems 
quenched across a critical point |3J. Examples include super- 
fluids |4'L, superconductors jstl, spin chains |2], cosmological 
models [6], and others [3]. 

In conventional superconductors, observations of the dy- 
namics of superconductivity would require a real-time change 
of either dopant concentration or lattice structure, injection 
of quasi-particles, or exposition to microwave or light 01, 
which is very difficult. The dynamics of BCS wavefunction 
is also of interest to nuclear physicists l^t]. Such dynamics 
can be realized in cold atoms by adjusting the inter-atomic 
interactions, which is routinely done in current experiments. 
Moreover, the dynamics of a BCS superfluid can be mapped 
to a central-spin problem 1 10]. Our goal is to study the non- 
equilibrium dynamics of an s-wave BCS superfluid resulting 
from a ramp of attractive pairing interactions. 

Previous work on the quench dynamics of a BCS super- 
fluid focused on the non-equilibrium gap followirm an instan- 
taneous change of the pairing interaction 111 ll4l3ll . Here we 
explore finite-rate quench dynamics which, unlike an instanta- 
neous one, allows for an observation of all stages of the evolu- 
tion ranging from adiabatic to strongly non-equilibrium ones. 
This can be done in a controllable way by changing the quench 
rate, which also allows for observations of universal critical 
properties of the phase transition through the scaling relations 
of out-of-equilibrium observables. These experimentally ac- 
cessible features should make a substantial progress in the un- 
derstanding of non-equilibrium physics of BCS superfluids. A 
linear ramp of the external magnetic field around the unitary 
^o^ t is predicted to drive the pairing gap out of equilibrium 

Several fundamentally important questions will be ad- 



dressed: What is the non-equilibrium gap and Cooper pair 
size during the ramp of the pairing interactions and how are 
the two quantities related? How do those non-equilibrium 
quantities reflect the critical properties characterizing the BCS 
superfluid-Fermi liquid (BCS-FL) phase transition? 

We consider a three-dimensional homogeneous two- 
component (cr I) unpolarized Fermi gas. Although in 
experiments there is a trap potential, the center of the Fermi 
gas can be considered as locally homogeneous and the dynam- 
ics we investigate should be observable. For weakly attractive 
or repulsive interactions the ground state is a BCS superfluid 
ifTsIl or a normal Fermi liquid |16], respectively. The Hamil- 
tonian with a time-dependent BCS-type pairing interaction is 

H = J2 '^kcljka - A(i) J2 cU^-ki^-il^n- (1) 

k(T kl 

Here = h'^k'^ /2m and X{t) is the time-dependent coupling 
constant from interactions between atoms in different hyper- 
fine states, which is regularized as A^^ = —m/Airh^a + 
^j,^^(l/2efc) iITtIi . a is the two-body s-wave scattering 
length and /C ^ A;f is a momentum cutoff ifisll . We choose 
as our units the Fermi momentum and energy kp and Ep 
of a non-interacting Fermi gas with the same density. Ex- 
perimentally a is tuned by an external magnetic field B via 
a = ai,g[l — 6B/{B — Bq)], where a^g is the background 
scattering length, while 6B and Bq are the width and position 
of the Feshbach resonance |15]. The T = BCS-FL quan- 
tum critical point corresponds to the point where the system 
becomes a non-interacting (free) Fermi gas. In BCS-Leggett 
theory the critical point is at a = and it can be located nu- 
merically for a finite-range potential 1 19]. This critical point 
for atoms is located at 209.96" IMQM- 

The excitation gap on the BCS side (o < 0) closes expo- 
nentially near the critical point. It is well approximated for 
kpa > -Iby iH 

A^i w 8£:Fe-''/2fcF|a|/g2^ ^2) 

Importantly, the excitation gap in BCS theory serves as the 
order parameter as well. This is a feature that distinguishes a 
BCS superfluid from other systems previously studied in the 
field of dynamics of quantum phase transitions f3\. 
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The BCS coherence length (after adjusting an overall con- 
stant) behaves near the quantum critical point as [TJ 



(3) 



Equations dU and (|3]l set the critical exponents to z = 1 
and v — oo through the relations ^ (A°'i)^^/^ and 



The infinite i' exponent makes a BCS super- 



fluid distinct from typical systems considered in the dynamics 
of phase transitions where ly — 0{l) |2]. To account for this 
singular behavior we will propose a quench protocol allowing 
for efficient tests of quench-induced scaling relations. 

To numerically study the effect of driving, we employ the 
generalized BCS-Leggett wavefunction 1 17] 



fct -fcJ. 



)|0>, 



(4) 



which provides a consistent description of BCS-Bose Einstein 
condensate (BEC) crossover but does not include medium ef- 
fects 1 11-13. 15. 17. 22], which renormalize the equilibrium 
gap. The medium effects can introduce corrections to the dy- 
namics, which may be evaluated non-perturbatively beyond 
BCS-Leggett theory [151 . Minimizing the quantum action 
/ dt{ih{^\dt\'if)/2 + c.c. - (^liJl^-)) iB, we obtain the 
following equations: 



hdk = 2igkA + 2idkek, hgk ^ idkA* 



idlA. 



(5) 



Here dk = 2ukv]:, gk — \vk\ ~ |ufe| , and the non- 
equilibrium gap function A (i) = A(i) (ifc(t)/2. Eqs. (|5]l 
are consistent with those in Ref. ifisll . Here tq = H/Ep is the 
unit of time. Eqs. (|4| and (|5]) can be appUed to the whole BCS- 
BEC crossover, though here we focus on the quench dynamics 
close to the a = critical point. 

The quench dynamics proceeds as follows. Initially the sys- 
tem is in its ground state at fci?a <C — 1, and the pairing inter- 
action is slowly decreased towards a = (the quantum critical 
point) by tuning away from the Feshbach resonance. In the be- 
ginning the system evolves adiabatically. At some point, say 
a — —a, the gap becomes so small that the system will fall out 
of equilibrium. This happens when the reaction time h/ A'^^ 
of the BCS superfluid becomes comparable to the inverse of 
the quench rate A^^i/lA^^I Q: 



A<'q(-a) 



|A<=q| 



(6) 



where C is a constant and the dot stands for d/dt. 

Basic understanding of what happens next comes from 
the Kibble-Zurek theory [3, 6], which has not been verified 
for a BCS superfluid yet. This theory predicts that a sys- 
tem switches suddenly from adiabatic to diabatic dynamics: 
its state is frozen near the critical point. Thus in the non- 
equilibrium stage of the evolution, the order parameter A 
and coherence length ^ should be given by their values A 
and ^ at the border between adiabatic and diabatic regimes: 
A = A'^'i(— a), ^ = a). In reality there is no sudden 
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Figure 1. (Color online) Non-equilibrium dynamics of the gap. Red 
solid line shows lA(t)| for tq — 150ro while black dashed line 
depicts A'^''. The two results start to deviate near kpait) — —0.75, 
which signals the onset of non-equilibrium dynamics. Upper inset: 
Illustration of how the threshold of non-equilibrium gap dynamics 
is determined. The solid green curve shows the relative difference 
between A'^'' and non-equilibrium gap for tq — 150to. When the 
relative difference reaches 50%, the time left to reach the critical 
point equals t. The power-law scaling of t is shown in the lower 
inset. The linear fit (solid red line) to numerical data (pluses) gives 
ln(£/ro) = -1.38±0.005-H(0.496±0.Q01) Htq/to)- It is in good 
agreement with the prediction i ~ tq^ i jlOi l. See fl8] for details. 



freeze out, but one can investigate if the following scahng re- 
lations hold near the critical point: 



A{t) ^ Af{t/i), m=i9it/i), 



(7) 



where / and g are unknown functions, t is the time left to 
reaching the critical point at the onset of non-equilibrium dy- 
namics, a{—t) = — a, and the critical point is reached at t = 0. 
Once the quench protocol, a{t), is specified, quantitative pre- 
dictions about system dynamics can be made. 

The linear ramp of the external magnetic field proposed in 
Ref. |fl4ll results in a non-linear drive of kpa{t) and makes it 
difficult to predict the dynamics. A more typical quench pro- 
tocol assumes a linear quench |2]. This coiTesponds here to 
having kpait) cx i/rg, where tq provides the quench time- 
scale: the slower the system is driven, the larger tq is. Such a 
quench can be realized by a non-linear (in time) sweep of Bit) 
near a Feshbach resonance. For a system with an exponen- 
tially closing gap, however, it introduces logarithmic correc- 
tions (in Tq) to the non-equilibrium length scale ^ |25], which 
demands very long quench times for verification of scaling re- 
lations. To overcome this problem, we propose the foUowing 
quench protocol 



A'\a{t))e^/^Ep = -t/rQ, 



(8) 



where time goes from ti < to such that kpaiti) < — 1, 
and the prefactor is for convenience. This can be inverted 
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for kpa > —1 using Eq. dU 

kFa{t) = ^/HtVT^)- (9) 

Such a ramp may be induced by B{t) = Bq + SB/[1 — 
Tr/kpUbg \n{t'^ /tq)]. To avoid inessential complications with 
the exact inversion of Eq. (O, we use the protocol shown in 
Eq. (|9]l through the whole evolution and choose the quench 
times so slow that the system is brought to kpa — —I nearly 
adiabatically. Therefore, significant excitations in our calcu- 
lations only happen at kpa > —1, where Eqs. (O and (|9]l 
become consistent. From Eqs. (|6]l and dH) one gets 



40 r 



kpa = 7r/ln(8CTQ/Toe ), A cx Jtq/tqEp 



kp^ oc \Jtq/to. (10) 

Eqs. ( fTOl i complement the scaling relations shown in Eqs. Q. 
A simple power-law dependence of A, t and ^ on the quench 
time scale will serve as a verification of our predictions. 

Typical dynamics of the gap function is illustrated in Fig. [T] 
Away from the critical point the gap is large enough to enforce 
adiabatic evolution in which A obtained from Eq. (|5]l matches 
/SP'^. As the gap becomes small enough, the relative deviation, 
(I A| — A'^'i) / A^"!, grows fast (see the upper inset of Fig. [T]!. 
From the simulations, this happens when the time left to reach 
the critical point, i, behaves as Tq^, which agrees with the 
prediction of Eq. ( fTOl i. Moreover, the difference | A| — A°'^ is 
controlled by the quench rate and increases as tq decreases. 
Finally, we have checked that the scaling relation for A near 
the critical point, Eq. (|7]i, holds. Very close to the a = 
critical point both A and /SP'^ approach zero. This reflects the 
vanishing coupling constant \{t) and does not imply that the 
system has reached an equilibrium state. 

Equilibrium gap can be inferred from radio-frequency 
(RF) spectroscopy. Modification of this technique to non- 
equilibrium systems has been proposed |26]. Therefore we 
expect that the non-equilibrium dynamics of the gap function 
described above will be measured experimentally. 

Next we study the Cooper pair size, which is defined as 



where 1^ stands for the Cooper pair wave-function: {k\4>) = 
Uk^k lll^l 12711 . In equilibrium the Cooper pair size ^^'^ is 
the same as the coherence length: ^""^ = S^'^'^. In a non- 
equilibrium situation it is interesting to address if the Cooper 
pair size encodes a non-equilibrium length scale and remains 
finite despite the fact that diverges near the critical point. 

Typical evolution of is depicted in Fig. |2] It con- 
sists of three consecutive stages. In the first one, the sys- 
tem is away from the critical point and it evolves adiabati- 
cally: ~ £.c'^- In the second one, the system enters non- 
equilibrium dynamics and so 7^ Cc'^- Interestingly, we find 
kp£,c ^ Ep/\A\ there. Thus, the relation between the non- 
equilibrium Cooper pair size and | A| is the same as that be- 
tween ^Jr"! and A°i. Importantly, it implies that the measure- 
ment of the non-equilibrium gap proposed in |26] can be im- 
plemented to estimate the non-equilibrium Cooper pair size. 
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Figure 2. (Color online) Dynamics of the Cooper pair size in the 
non-equilibrium quench with tq = 150ro. The curves, from left 
to right, correspond to the non-equilibrium Cooper pair size kp^c 
(black), the inverse of the nonequilibrium gap, oc Ef/\A\ (red), and 
the instantaneous equilibrium pair size kp^c'^ (green). The blue box 
highlights the region where 7^ £,c'^ but oc Ef/\A\. The pro- 
portionality factor for the middle (red) curve was chosen to match 
the static result kpa = —4 (not shown on the plot). See fisll for 
details. 



In the third stage both kp^^'^ and Ep/\A\ diverge, but the 
non-equilibrium pair size stays finite. Now decouples from 
the gap function and depends solely on the non-equilibrium 
scales ^ and t through the relation shown in Eq. (|2|, which is 
depicted in Fig. [3] The inset of Fig. [3] quantitatively confirms 
that the non-equilibrium length scale ^ persists in the system 
till the end of evolution and scales as Tq^, which agrees with 
the prediction of Eqs. ( fTOb . Surprisingly, the existence of the 
two stages of non-equilibrium dynamics is not predicted by 
the Kibble-Zurek theory. 

How to measure the non-equilibrium Cooper pair size is an 
important open question. There are experiments "measuring" 
Cooper pair size by inverting the gap obtained from RF spec- 
troscopy ll22ll . This will not work close to the critical point 
(third stage of evolution on Fig.|2]) where and A are decou- 
pled. One possibility is to generalize to a non-equilibrium set- 
ting a scheme from 1I28I1 . where is inferred from the damp- 
ing rate of the collective mode of a coexisting BEC. 

We have showed that close to the BCS-FL quantum phase 
transition at T = 0, BCS superfluids exhibit multi-stage dy- 
namics with scaling behavior consistent with Kibble-Zurek 
theory. At finite temperature two additional effects show up. 
First, there is a classical (thermal) phase transition from a BCS 
superfluid to a normal Fermi gas at ksTc ~ Q.57A°'^{T — 0) 
[7], where fc^ is the Boltzmann constant. Second, quasi- 
particles (QPs) from broken Cooper pairs exist at T > 0. 
Their relaxation time is estimated to be Tqp « hEp / {A'q-^}'^ 



1 1 in . where the Ginzburg -Landau gap function Aq^ corre 



sponds to the effective gap in the presence of both QPs and 
Cooper pairs. Below we assume that the temperature changes 
insignificantly during the slow driving and estimate the influ- 
ence of these finite temperature effects. 

If we start the quench dynamics from an equilibrium state 
at T < Tc(a) and /cfa < 0, Aq\ « A'^i initiaUy and the 
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Figure 3. (Color online) Scaling of the non-equilibrium Cooper 
pair size, ^c, for tq = 150ro (black line) and for tq = 3360to 
(red dots). Both curves collapse into a universal one near the quan- 
tum critical point, which confirms that the quench-imprinted non- 
equilibrium length scale ^ has not been "erased" until the end of the 
evolution (t — >■ 0). Rescaling was done with kpi, = \/ tq /to and 
t = y^TQTo. Inset: at the critical point, a{t = 0) = 0. Numerical 
data are shown as pluses. They are fitted by ln(fcF^c(a = 0)) = 
0.861 ± 0.002 + (0.4995 ± 0.0003) ln(rQ/ro) (green line). The 
expected scaling is ^c{a = 0) ~ tq^. See H8i1 for details. 



quasi-particle relaxation time Tqp will be short compared to 
the inverse of the quench rate A^^/I^glI- Thus both the su- 
perfluid and QPs follow the ramp adiabatically. Comparing 
the two time scales as we did in Eq. QPs will stop fol- 
lowing the quench dynamics of Eq. ^ when the time left to 

reach the quantum critical point is t* ^ ^TqTq. Comparing 
this to the onset of the non-equilibrium dynamics of the con- 
densate, t ^ ^tqTq, one can see that for tq/tq ^ 1 studied 
in this paper the QPs will fall out of equilibrium before the su- 
perfluid does {t* > t). Here we assume ~ ^""^^ which 
should be reasonable if the system is away from the classical 
critical point at time t* . This requires that T <C Tc{t*) ^ 
^tq/tqTf, where Tp = Ep/kB- After QPs fall out of 
equilibrium, they approximately decouple from the superfluid 
dynamics. Next the superfluid falls out of equilibrium. We 
expect that the influence of the quantum critical point on the 
dynamics (e.g. through the Cooper pair size dependence on ^) 
will be visible when T ^ Tf.{i) ^ \/tq/tqTp. This ensures 
that there is still a condensate at the time when quantum criti- 
cality drives the superfluid out of equilibrium. Experimentally 
one can stop the ramp before the system reaches the estimated 
classical transition line and observe the predicted scaling from 
a quantum quench. 

Since the equilibrium Tc vanishes as a 0, the system 
will eventually be swept out of the BCS superfluid phase. 
This raises another important question: Is the BCS superfluid- 
Fermi gas classical phase transition still well defined when 
both QPs and the superfluid are out of equilibrium? We point 
out that crossing a classical critical point may also imprint 



additional non-equilibrium length scale encoding the critical 
exponents of the classical transition [3]. 

In summary, we have shown that a finite-rate quench al- 
lows for controlled (through the quench time scale tq) stud- 
ies of non-equilibrium dynamics of a BCS superfluid. The 
quench imprints non-equilibrium energy and length scales vis- 
ible through the gap and Cooper pair size. Our results call for 
further development of experimental techniques for studies of 
non-equilibrium cold Fermi gases and point out that finite- 
temperature effects may induce additional scaling behavior of 
non-equilibrium BCS superfluids in a finite-rate quench. 

This work is supported by U.S. Department of Energy 
through the LANL/LDRD Program. 
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